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Abstract. Let o~t(x) denote the implied volatility at maturity t for a strike K = Soe , where 
i£t and So is the current value of the underlying. We show that o~t(x) has a uniform (in a;) limit 
as maturity t tends to infinity, given by the formula a oc (x) — y/2 ^/t*(x) 1//2 + (h*(x) — x) 1 ^ 2 ^ , for 
x in some compact neighbourhood of zero in the class of afhne stochastic volatility models. The 
function h* is the convex dual of the limiting cumulant generating function h of the scaled log-spot 
process. We express h in terms of the functional characteristics of the underlying model. The proof 
of the limiting formula rests on the large deviation behaviour of the scaled log-spot process as time 
tends to infinity. We apply our results to obtain the limiting smile for several classes of stochastic 
(*V^ volatility models with jumps used in applications (e.g. Heston with state-independent jumps, Bates 

Qh ' with state-dependent jumps and Barndorff-Nielsen-Shephard model). 

d " 
<+- 1 

1. Introduction 

Let the process S = e model a risky security under an equivalent martingale measure and let 

OO ! 

at(x) denote the implied volatility at maturity t for a strike K = Soe (see (j50jl for the precise 



definition of at(x)). The main result of the present paper (Theorem |14p states that, if the log-spot 
X follows an affine stochastic volatility process with jumps, then o*t(x) converges to <r 00 (x) as the 



maturity t tends to infinity, where <7oo is given by the formula 

(1) ^oo(x) = ^{h*^) 1 ' 2 + {h*(x)- X ) 1/2 ). 

• r- I . 

The function h is the limiting cumulant generating function of the scaled log-spot (Xt/t)t>\ and 
h* is its convex dual (i.e. the Fenchel-Legendre transform of h). Locally uniform convergence of 
the implied volatility to (7^ is also established. 

In [FJlll [FJM11] the limiting behaviour of the smile at large maturities in the Heston model is 
investigated. Theorem 1141 can be viewed as a generalisation of the main result in |FJ1H IFJM U]. 
Not only does it cover a large class of stochastic volatility models with jumps rather than a single 
affine model with continuous trajectories, but furthermore provides a better understanding of the 
limit: Theorem [TH states that the limit holds also at the critical points x* and x*, which are 
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excluded from the analysis in |FJ11[ fFJMll]. and the convergence on the set R\ {x*,x*} is shown 
to be uniform on compact subsets. 

In the class of affine stochastic volatility models, the formula for the limiting implied volatility 
for a fixed strike proved in Tehranchi [Teh09] (see also [LewOO] in the case of the Heston model) also 
follows from (pQ), since in Theoremll4lthe convergence is uniform on a compact neighbourhood of the 
origin. In [GLllj . the authors give various representations for the implied volatility, including in the 
large-maturity regime, based on an assumed asymptotic behaviour of certain European derivatives 
in the underlying model, which is not specified. This representation is not fully explicit in terms 
of the model parameters and it is therefore unclear how to apply it directly to the class of affine 
stochastic volatility models. 

Contribution of the paper is twofold. First we study the properties of the limiting cumulant 
generating function h of the affine stochastic volatility models. Results in Lemma [9j Theorem [TU1 
and Corollary [TT] give new properties of the function h, which are crucial for the understanding 
of the large deviation behaviour of the model. Second, the problem of understanding the limiting 
behaviour of option prices and the corresponding implied volatilities using the large deviation 
principle is tackled. The uniform limit in x (on all compact subsets of M) of vanilla option prices 
is given in Theorem 1131 for non-degenerate affine stochastic volatility models and exponential Levy 
models (i.e. degenerate affine stochastic volatility models). As mentioned above Theorem [T4l deals 
with the limiting implied volatility smiles in these classes of models. 

Besides giving a formula, which relates model parameters and the limiting implied volatility 
smile, these theoretical results have the following practical consequences: 

(1) in the large-maturity regime studied in this paper, the jumps in the model influence the limiting 
implied volatility smile as maturity tends to infinity (see examples in 16. ID : 

(2) for every affine stochastic volatility model there exists an exponential Levy model such that 
the smiles of the two models in the limit coincide. In other words the stochasticity of volatility 
does not (in the affine class) enlarge the family of possible limiting implied volatility smiles (see 
Section [U] for details). 

The starting point of the analysis of the large deviation behaviour of an affine stochastic volatil- 
ity process (X, V) in the present paper is Theorem [HI taken from jKRlll Theorem 3.4]. This result 
describes certain properties of the limiting cumulant generating function h, which are however in- 
sufficient to understand the essential smoothness of h required in establishing the large deviation 
principle of (Xt/t)t>\- The main contribution of this paper in the area of affine processes is Theo- 
rem [TU] which identifies sufficient conditions for the process (X, V) that imply essential smoothness 
of the function h. The conditions in Theorem [10] are easy to apply to the models of interest (see 
e.g. Section [22]) . Its proof goes beyond the analysis in |KR11| as one is forced to study the special 
Levy-Khintchine form of the characteristics of the process, since their general convexity properties 
no longer suffice to establish the required behaviour of the limit. 
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The rest of the paper is organized as follows. In Section 2 we define the class of affine stochastic 
volatility processes and recall some of their properties. In Section 3 we review briefly basic concepts 
in the theory of large deviations and state the Gartner-Ellis theorem. Section 4 establishes the large 
deviations principle for the scaled log-stock of an affine stochastic volatility model as maturity tends 
to infinity. Sections 5 and 6 respectively translate this result into option price and implied volatility 
asymptotics. Numerical examples are given at the end of Section 6. 

2. Affine stochastic volatility models with jumps 
Consider a stochastic model for a risky security S = (St)t>o given by 
(2) S t = exp((r - d)t + X t ), t>0, 

where the interest rate r and the dividend yield d are non-negative and constant and the log-price 
process X = (X t )t>o starts at Xq € R. Since the dynamics of S is given under a risk- neutral 
measure, the forward price process is (exjp(X t ))t>o- We assume throughout the paper without loss 
of generality that S is a forward price process (i.e. r = d). Denote by V = (Vj)t>o a process, starting 
at a constant level Vo > 0. The process V can be interpreted as the instantaneous variance process 
of X but may also control the arrival rate of jumps of X. We make the following assumptions on 
the process (X, V) throughout the paper. 

Al: (X, V) is a stochastically continuous, time-homogeneous Markov process with state-space 

D = R x M >0 , where R^ := [0,oo). 
A2: The cumulant generating function &t( u , w) of (Xt, Vt) is of a particular affine form: there 

exist functions 4>(t, u, w) and ip(t, u, w) such that 

$ t (u, w) := log E [exp(uX t + wV t )\ X , V ] 
= <f>(t, u, w) + Voip(t,u, w) + X u 

for all (t,u,w) S M^o x C 2 , where the expectation exists. 

Remarks, (i) Al and A2 make (X, V) into an affine process in the sense of [DFS03 . 

(ii) Al and A2 imply a homogeneity property of (X, V): if the starting value (Xq, Vq) is shifted by 
(x, 0) G D, the law of the random variable (Xt, Vt) is shifted by the vector (x, 0) for any t > 0. 

(iii) Assumptions Al and A2 imply that the variance process V is a one-dimensional strong Markov 
process in its own right. 

(iv) The law of iterated expectations applied to $t(u, w) yields the flow-equations for 4> and ij) (see 
[DFS031 Eq. (3.8)-(3.9)]): 

4>(t + s, u, w) = (f)(t, u, w) + (f>(s, u, ip(t, u, w)), 

(3) 

ip(t + s, u, w) = tp(s, u, tp(t, u, w)), 

for all t, s > 0. 
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(v) It is shown in |KB.1U Thm. 2.1] (see also [DFS03]) that if \<f)(T,u,ri)\,\tp(T,u,rf)\ < oo for 
(t,u,t]) G (0, oo) x C 2 , then for all t G [0,r) and w £ C such that 5Rw < Sir/, the functions <f> and 
ip satisfy the generalized Riccati equations 

(4a) dt4>(t,u,w) = F(u,ifi(t,u,w)), </>(0, u, w) = 0, 

(4b) d t if)(t,u,w) = R(u,ip(t,u,w)), ip(0,u,w) = w, 



where 

d 

(5) F(u,w) := —4>(t,u,w) 



d 

R(u,w) := —ip(t,u,w) 
t=o+ ot 



t=o+ 



Furthermore for all t G [0,r] we have \(f)(t,u,w)\, \ip(t,u,w)\ < oo. 

(vi) If (X, V) is a diffusion process, then ODEs fl3J) become classical Riccati. Note also that 
follows from the flow equations ([3]) by differentiation with respect to s. 

(vii) (j> and ^ can for small t be expressed implicitly in terms of F and R as 

/•* ril>(t,u,v>) ^„ 

d)(t,u,w) = / F(u,ib(s,u,w)) ds and / — r=i- 

Jo 7i« 

The functions i 7 and i?, defined in ([5]), must be of Levy-Khintchine form (see [DFS03 ). In other 
words 

(6a) F(u, w) = (f 0", ™)'> C"> ^ )') + ( fc > 0> W Y) ~ c 

+ [ (e<^)'> - 1 - (lj f (0, (u, w)')) m(d£), 
Jd\{o} v 7 

(6b) w) = (|(u, to)', (u, «;)') + (P, (u, «;)') - 7 

+ / (e<^ w ^'> - 1 - (cj r (0, (u,w)')) M (d0, 

where D = lx R>0j (v) is the inner product on R 2 , (u,w)' denotes transposition, up, ujr are 
suitable truncation functions, which we fix by defining 

MO = and = ^ j , where £ = (^j , 

and the parameters (a, a, 6, (3, m, n) satisfy the following admissibility conditions: 

• a, a are positive semi-definite 2 x 2-matrices with ai2 = <22i = ^22 = 0; 

• 6 G D, /3 G M. 2 and c, 7 G R >0 ; 

• m and \i are Levy measures on D and /o\{ } ((£1 + £2) A l) m(d£) < 00. 
Assumptions Al and A2, the generalized Riccati equations and the Levy-Khintchine decompo- 
sition ([6]) lead to the following interpretation of F and R: F characterizes the state-independent 
dynamics of the process (A, V) while R characterizes its state-dependent dynamics. The instanta- 
neous characteristics of the Markov process (X, V) are given as follows: a + Va the instantaneous 
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covariance matrix, b + V (3 the instantaneous drift, m(d£) + V^i(d^) the instantaneous arrival rate 
of jumps with jump heights in d£ and c + jV the instantaneous killing rate. 

The function \ defined below plays a key role in the characterisation of the martingale property 
of the process S = exp(X). 

Definition 1. For each m£R such that R(u,0) < oo, define x( u ) as 



Remarks, (i) The condition R(u,0) < oo implies that, for some 5 > the function w \— > R(u,w) is 
convex on (—5,0] and differentiable on (—5,0), since the process V does not have negative jumps. 
Therefore x( u ) is a well-defined, possibly equal to +oo, convex function given by the limit of 
&2R(u, w) as w t 0. It can be expressed explicitly as 



(ii) The sufficient and necessary condition for S to be conservative and a martingale, in terms of 
R,F and Xi is given in |KR1H Thm. 2.5]. A simple sufficient condition for these properties reads 
(see [KRlll Cor. 2.7]): 

• if 2^(0, 0) = 22(0, 0) = and x(0) < oo then S = exp(X) is conservative. 

• if 5* is conservative, F(1,Q) = 22(1,0) = and x(l) < oo, then S = exp(X) is a martingale. 

Since S serves as a forward price process under a risk-neutral measure P in an arbitrage-free 
asset pricing model, it has to be conservative and a martingale and hence we assume: 

A3: F(0, 0) = R(0, 0) = F(l, 0) = 22(1, 0) = and x(0) + x(l) < oo. 
In particular 2^(0, 0) = 2?(0, 0) = in assumption A3 means that the instantaneous killing rates 
c and 7 in © are zero and the condition F(1,0) = 22(1,0) = is closely related to the functions 
ip(-, 1,0) and 4>(-, 1,0) being identically equal to zero (see the generalized Riccati equations in Q), 
which implies the martingale property of S = exp(X). The following non-degeneracy assumption 
will guarantee the stochasticity of volatility of the process X. 

A4: There exists some u E M, such that R(u, 0) 7^ 0. 

Definition 2. The process (X, V) is a non- degenerate (resp. degenerate) affine stochastic volatility 
process if it satisfies assumptions Al - A4 (resp. Al - A3 and does not satisfy A4) and S = e x is 
the corresponding affine stochastic volatility model. 

Remark. Assumption A4 excludes the degenerate case where the distribution of X does not de- 
pend on the volatility state Vq. Indeed, if A4 is not satisfied, i.e. 2?(-,0) = 0, then ^ implies 
that ip(t, u,0) =0 and <j)(t, u, 0) = tF(u, 0) for all (t, u) E M^o x C where the expectation in as- 
sumption A2 exists. Hence if A4 does not hold, then A2, (|6ap and the characterisation theorem for 
regular affine processes |DFSQ3[ Theorem 2.7] imply that S = e x is an exponential Levy model. 



X{u) := d 2 R{u,w)\ 



w=0 ■ 




w=0 




where £' = (fi,6)- 
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In particular the class of affine stochastic volatility models includes the Black-Scholes model as a 
degenerate case. 

The following proposition describes certain properties of F and R that will play a crucial role in 
Section gTj 

Proposition 3. Let (X, V) be a non- degenerate affine stochastic volatility model and let the sets 
T>f = {( u , w) €zM? : F(u, w) < oo } and V R = {(u,w) G M 2 : R(u, w) < 00} be the effective do- 
mains of the functions F and R respectively. Then the following holds: 

(A) F and R are lower semicontinuous convex functions, which are continuously differentiate in 
the interiors T>° F and T>° R (in M. 2 ), and their effective domains Dp and T>r are also convex; 

(B) F and R are either affine or strictly convex functions when restricted to one- dimensional affine 
subspaces ofM 2 . 

Proof. The Levy-Khintchine representation for F and R in © implies that they are cumulant 
generating functions of some (infinitely divisible) random vectors taking values in M 2 . Holder's 
inequality yields that F and R are convex. The dominated convergence theorem and the represen- 
tation in ([6]) implies that F and R are analytic in Dp and T> R respectively. Fatou's lemma implies 
that the functions F and R are lower semicontinuous. Since F and R are cumulant generating 
functions, the second derivative of their restriction to an affine subspace is either identically zero 
or strictly positive everywhere (each affine subspace in M? corresponds to a random variable which 
takes values in M and may or may not be constant almost surely). This concludes the proof. □ 

2.1. SDE representation of affine stochastic volatility processes. In order to define an affine 
stochastic volatility model one needs to choose admissible parameters (a, a, b, (3, m, fi) such that the 
corresponding process (X, V), which exists by |DFS031 Thm. 2.7], satisfies assumptions Al - A3 
(note that c = 7 = by A3 and will henceforth be ignored). This procedure yields a semigroup, 
and hence the law, of the Markov process (X, V) which is in principle sufficient for option pricing. 
However path-wise descriptions of the pricing models in financial markets are widely used as they 
add to the intuitive understanding of the properties of the model. In the rest of this section we 
briefly describe a path- wise construction of the process (X, V), given in |DL06] . and relate it to the 
most popular affine stochastic volatility models used in derivatives pricing. 

Assume that the parameters (a, a, b, f3, m, ^1) are admissible and suppose in addition that the 
tails (i.e. the large jumps) of m and \x satisfy: 



(7) 



/ ||£||m(d£)<oo and / ||£|| M (d£) < 00 
J\\i\\>^ hz\>t 
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where ||f|| = (£,£}. Let 



(8) 



bi 



L 
L 




^m(d£), b 2 = b 2 , 



~2 A*(d£) for i = 1,2. 



D\{0} 1 + £j 



Note that the integrals in (JSj) are finite by (|7|) and the parameters (a, a, 6, /3, m, fx) are admissible 
with appropriate truncation functions. 

Let (J r t)t>o, P) be a filtered probability space equipped with 

• a three-dimensional standard Brownian motion (B , i? 1 , i? 2 ), 

• a Poisson random measure A^o(ds,d^) on IR^o x D with compensator dsm(d^), 

• a Poisson random measure Nx(ds, du, d£) on R> x D with compensator dsdu/i(d^), 

where as usual D = R x Mj, denotes the state-space of the model. Let 

iV (ds, d£) = iV (ds, df ) - ds m(d£) and JVi(ds, du, d£) = iV~i (ds, du, df ) - ds du u(d£) 

be the compensated Poisson random measures and let a be a 2 x 2-matrix such that aa T = a. 
Theorem 6.2 in [DL06J implies that the system of SDEs 



with initial condition (Xq, Vq) £ M x (0, oo), has a unique strong solution (X, V) that is an affine 
Markov process with admissible parameters (a, a, b, (3, to, u). 

Remarks, (i) The change of parameters b and (3 introduced in (JSj) is inessential. Its function is to 
establish the notational compatibility with [DL06J. 

(ii) The integrals in ([9jH (|10p against N\ are taken over a random set whose dsdu«(d£)-volume is 
proportional to This, together with the structure of the Poisson random measure N\, reinforces 
the intuition that the jumps of the process (X, V) that correspond to the integral term with respect 
to N\ have random intensity which is proportional to V. 

2.2. Examples of affine stochastic volatility models. We now describe some of the affine 
stochastic volatility models that are of interest in the financial markets and can be obtained as 
solutions of the special cases of SDE (f9^- (fT0j) , 



(10) 



(9) 
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2.2.1. Heston model. The log-price X and the stochastic variance process V are given under the 
risk-neutral measure by the SDE 

dX t = -^dt+y/VtdW?, 

dV t = -X(V t -9)dt + Cv^dWf , 

where W l ,W 2 are Brownian motions with correlation parameter p E (—1,1), and C>^>0 > 
(see [Hcs93j). The affine characteristics of the model are 



(11a) F(u, w) = X6w, 

1 c 2 

(lib) R(u, w) = — (u 2 — u) + -^w 2 — Xw + uwpC- 

It is easily seen that \ is given by 

(12) X (u) = p(u - X 



and it is trivial to check that Al - A4 are satisfied. 

2.2.2. Heston model with state-independent jumps. Let J be a pure-jump Levy process independent 
of the correlated Brownian motions W 1 and W 2 . The Heston- with-jumps model is defined by the 
SDEs 

dx t = (s-^]dt+ Vvidwl + dJ u 

dV t = -X(V t -0)dt + cVVtdw?, 

where £, A, 9 > and 5 G M. Assume that J is a spectrally negative Levy process with characteristic 
exponent jlogE[e uJt ] = /(.^o^e^ 1 — 1 — u£i/(l + ^f))^(d^i). Since J only jumps down, this 
captures the generic situation in the modelling of equity markets. Assume further that the jumps 
of J are integrable (i.e. C^o) l^ 1 ! u (^ 1 ) < °°)- ^ n or der to identify the coefficients in (j9])- (fT0|) . 
we compensate J so that it becomes a martingale and can hence be expressed as an integral against 
the compensated Poisson random measure Nq. This implies 

bi = 6+[ T^Kdei). 

J (-oofl) 1 +4l 

It is easily seen that a = 0, an = 1, ct\2 = pC, «22 = C 2 i M = and m(d£) = [y ® #o)(d£), where 
<5q is the Dirac delta measure. Therefore b% = 5, 62 = /?i = —1/2 and fa = —X. The martingale 
condition (F(1,0) = 0) implies 5 = — fr^ \ (e^ 1 — 1 — £i/(l + £f)) f(d£i) and the affine form of 
the model is given by 

(13a) F(u,w) = X6w + k(u), 

1 c 2 

(13b) R(u, w) = -(u 2 — u) + —w 2 — Xw + uwpC, 
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where k(u) is the compensated cumulant generating function of the jump part, i.e. 

(14) k{u) = [ fe fin - 1 -u(e^ - l)) i/(dft). 

J (-oo,0) v \ J J 

2.2.3. A model of Bates with state- dependent jumps. We consider the model given by 
dX t = -(l + s]v t dt+ VVtdW? + f Ci N(V t , dt, d&), 

dT4 = -A(Vt - 0) dt + C V^dW?, 

where as before A, 6, £ > 0, 5 G R and the Brownian motions VF 1 and VF 2 are correlated with corre- 
lation p 6 (—1,1). The jump component is given by JV(Vt, dt, d£i) = N(Vt , dt, d£i) — n(T4, dt, d£i), 
where A^Vj, dt, d£i) is a Poisson random measure independent of W 1 and W 2 with intensity mea- 
sure n(Vt, dt, d£i) of the state- dependent form V^(d£i)dt. Here f(d£i) denotes a Levy measure 
onI\ {0}. A model of this kind has been proposed in [BatOOj to explain the time- variation of 
jump-risk implicit in observed option prices. 

As in Section 12.2.21 we assume that the support of v(d£i) is contained in (— oo,0) and that 
the inequality C^q) l^ 1 ! K^Ci) < oo is satisfied. We can identify the parameters in (f9|)- (fT0]) as 
a = 0, an = 1, ai2 = p(, a 22 = ( 2 , Pi = -1/2 - 8, /3 2 = -A, b\ = 0, b 2 = b 2 = \0, m = and 
p{d^) = (y ® 5o)(d^), where So is the Dirac delta concentrated at 0. Hence we find 

ft = 4 -S- f TT^^ 

1 J (-oo,0) 1 +?i 



and /?2 = — A, &i = 0. The functions -F and i? for the Bates model are 
(15a) F(u,w)=\6w, 

1 c 2 

(15b) R(u, w) = — (u 2 — u) + -^w 2 ~ ^ w + uwpC + k(u), 

where k(u) = f/^ \ (e^ lU — 1 — u (e^ 1 — l)) f (d£i) and the martingale property (R(l, 0) = 0) was 
used to determine the value of the parameter S = // ^ ) ~~ ^ ~~ ^0 K^Cl)- ^ is clear that 
x(^) = pC'W ~~ A and that Al - A4 are satisfied. 

2.2.4. TTte Barndorff-Nielsen-Shephard (BNS) model. The BNS model was introduced in |BNS01j 
as a model for asset pricing. Under a risk-neutral measure, it can be defined by the following SDE 

dX t = (5- ±V t )dt + y/V t dW t + pdJ xt , 
dV t = -XV t dt + dJ xu 

where A > 0, p < and (Jt)t>o is a Levy subordinator with the Levy measure v, i.e. a pure jump 
Levy process that increases a.s. The cumulant generating function k(u) of (Jt)t>o takes the form 

(16) «(«) = ( (e^ 2 - 1) i/(d&). 

■/(0,DO) 
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To conform with (J7]) we further assume that J*/ Q s £2 ^(d^) < oo. The drift 5 will be determined 
by the martingale condition for S. The time-scaling J\t is introduced in [BNS01J to make the 
invariant distribution of the variance process independent of A. The distinctive features of the BNS 
model are that the variance process has no diffusion component, i.e. moves purely by jumps, and 
the negative correlation between variance and price movements is achieved by simultaneous jumps 
in V and X. 

It follows from ([§|)- (|10p and the SDE above that a = ct\2 = 022 = 0, an = 1, p = and 

m(dO = I { £ 1= p 6} Ai/(d£ 2 ), 

where I^ 1=p ^ denotes the indicator function of the half- line £1 = p& in Z? \ {0}. Therefore it 
follows that p 1= p 1 = -1/2, p 2 = fo = -A, b 2 = b 2 = 0, h = S + A J (0 oo) p£ 2 i>(d&) and 

J (0,00) 1 + (p&r 

The definition of F in ([6]) and the martingale condition F(1,0) = imply that we need to define 
5 = —Xk(p), where k is the cumulant generating function of J given in (|16p . The BNS model is an 
affine stochastic volatility model with F and R given by 

(17a) F(u, w) = Xk(w + pu) — u\n(p), 

(17b) R(u,w) = i(u 2 -u) - Xw. 

We have x( u ) = — A and the assumptions Al - A4 are clearly satisfied. 



3. Large deviation principle and the Gartner-Ellis theorem 

In this section we give a brief review of the key concepts of large deviations for a family of 
(possibly dependent) random variables (Zt)t>i and state a version of the Gartner-Ellis theorem 
(see Theorem ED that will be used to obtain the asymptotic behaviour of the option prices and 
implied volatilities. A general reference for all the concepts in this section is [DZ98, Section 2.3]. 

Let Zt take values in R and recall that I : R — > (—00, 00] is lower semicontinuous if {x : I (x) < a} 
is closed in R for any a G R (intuitively for any xq E R the values of I near xq are either close 
to I(xq) or greater than I(xq)). A nonnegative lower semicontinuous function / is called a rate 
function. If in addition {x : I (x) < a} is compact for any a G R, then / is a good rate function. 

Definition 4. The family (Zt)t>i satisfies the large deviation principle (LDP) with the rate func- 
tion I if for every Borel set BcRwe have 

- inf{/(x) : x £ B°} < liminf - log P [Z t EB]< limsup - log P [Z t G B] < - inf \l(x) : x G B~\ , 

£-►00 t t^-oo t 

with the convention inf = 00 (the interior B° and closure B are relative to the topology of R). 
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An important consequence of Definition is that if (Zt)t>i satisfies LDP and I is continuous on 
B, then lim^oo t~ l log P [Z t G B] = - inf{/(x) : x G B}. 

The Gartner-Ellis theorem (Theorem[6]) gives sufficient conditions for {Zt)t>\ to satisfy the LDP 



where T>\ :={«£! : A(u) < 00} is the effective domain of A and T>°^ is its interior in R. Since Af 
is convex (by the Holder inequality) for every t, the limit A is also convex by |Roc701 Theorem 10.8] 
and the set T>\ is an interval. Since A(0) = 0, convexity of A and G T>°^ imply A(u) > —00 for all 
u G R. Furthermore the convexity implies that A is continuous on T>\. The statement in (|18p is 
an important assumption of Gartner-Ellis theorem (Theorem [6] below), which in particular implies 
2? A 7^ However the converse does not hold in general, i.e. if is a boundary point of a domain 
T>\ with non-empty interior, LDP may still hold true. 

A further property of the function A : R — > (—00, 00], which arises as an assumption in TheoremEl 
is essential smoothness. 

Definition 5. A convex function A : R — > (—00,00] is essentially smooth if 

(a) 2? A is non-empty; 

(b) A is differentiate in V° K \ 

(c) A is steep, in other words it satisfies lim n _ ) . 0O |A'(n n )| = 00 for every sequence (u n ) rag N in 
T>\ that converges to a boundary point of V\. 

The Fenchel-Legendre transform (or convex dual) A* of A is defined by the formula 



with an effective domain T>\* := {x G R : A*(x) < 00}. The following properties are immediate 
from the definition: 

(i) < A*(x) < 00 for all x G R, since A(0) = 0; 

(ii) A*(x) = supjiix — A(u) : u G T>\} for all x £ R and hence A* is convex in the interval Pa* 
and continuous in the interior "D A „ ; 

(iii) A* is lower semicontinuous on R as it is a supremum of continuous (in fact linear) functions. 
Hence the level sets {x : A*(x) < a} are closed. 

In general T>\* can be strictly contained in R and A* can be discontinuous at the boundary of T>\* 
(see |DZ98t Section 2.3] for elementary examples of such rate functions). Assumption (118p implies 
that for any 5 > 0, such that (—5, 5) C V\, and c = sup{A(n) : u £ [—5, 5]} we have 



Hence the set {x : A*(x) < a} is compact for any «£l and therefore A* is a good rate function. 




(19) 



A*(x) 



supjnx — A(u) : u £ R} for x G R 



(20) 



A*(x) > sup{ux — A(u) : u G [— 5, 5]} > 5\x\ — c. 
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Remarks. (A) If A is strictly convex, differentiable on D A and steep, which is the case in the 
applications in this paper, then D\* = R and for each x G R the equation A' (it) = x has a unique 
solution u x in T>°^. Furthermore the formula 

(21) A*(x) = xu x -A(u x ) 

holds. This reduces the computation of A*(x) to finding the unique root of the equation A' (it) = x, 
where the strictly increasing function A' is in most applications known in closed form. 
(B) If (Zt)t>i satisfies f)18|) and the function A satisfies the assumptions of Remark (A) and is twice 
differentiable with A" (it) > for all it G T>^, then ([2ip implies that the Fenchel-Legendre transform 
A* is differentiable with the derivative 

(22) (A*)'(x) = (A) -1 (a;) for all x G R. 

In particular (|22p implies that A* is strictly convex on R and that its global minimum is attained 
at the unique point x* given by 

x* = A(0). 

We state a simple version of the Gartner-Ellis theorem (for the proof see [DZ98j Section 2.3]). 

Theorem 6. Let (Zt)t>i be a family of random variables that satisfies assumption (I18D with the 
limiting cumulant generating function A : R — > (—00,00]. If A is essentially smooth and lower 
semicontinuous, then the LDP holds for (Zt)t>i with the good rate function A*. 

4. Limiting cumulant generating function in affine stochastic volatility models 

4.1. Non-degenerate affine stochastic volatility processes. Let (X, V) be a non-degenerate 
affine stochastic volatility process (see Definition [2]) . The goal of the present section is to describe 
the limiting cumulant generating function h of the family of variables (Xt/t)t>i, defined by 

(23) h(u) := lim - log E \e uXt ] 

t— >oo t 

for every it G R where the limit in (|23[) exists as an extended real number. The function h will 
determine the limiting implied volatility smile of the model S = e . To ensure that h is finite on 
an interval that contains [0, 1] , which is key for establishing the LDP, a further assumption will be 
required: 

A5: x(0) < and x(l) < 0, where x is given in Definition [TJ 
This assumption will also imply that h can be uniquely extended to a cumulant generating function 
of an infinitely divisible random variable. 

In order to apply the Gartner-Ellis theorem in our setting, we need to answer the following three 
questions: is h well-defined as an extended real number by (f2"3"j) for every it G R, does the effective 
domain 2\ contain [0, 1] in its interior and is h essentially smooth? Answers to these questions 
play a crucial role in establishing the large deviation principle, via Theorem [61 for affine stochastic 
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volatility models. Theorem [10] and Corollary [TT| proved in this section, provide easy to check 
sufficient conditions for the affirmative answers to hold. 

It is shown in [KRllj that the function h can be obtained from the functions F and R without the 
explicit knowledge of (ft an d rp (see Section [2] for definition of Lemma [7] and Theorem [8] taken 

from [KRTTJ Lemma 3.2 and Theorem 3.4], describe certain properties of the limiting cumulant 
generating function h, which are needed in Section [5] but are insufficient to guarantee the essential 
smoothness of h. The main contribution of the present section is Theorem [Tu] which identifies 
sufficient conditions for the process (X, V) that imply essential smoothness of the function h. The 
conditions in Theorem [10] are easy to apply to the models of Section 12.21 which will allow us to find 
their limiting implied volatility smiles. 

Lemma 7. Let (X, V) be a non- degenerate affine stochastic volatility process that satisfies assump- 
tion A5. Then there exist a maximal interval X and a unique convex function w : X — > M such that 

w e C(i)nC(r) and 

R{u,w{u)) = for all u £ X, 

where R is given in ([5]) (see also ©j. Furthermore we have 

(a) [0, 1] C X and d 2 R(u, w(u)) < for all u £ X° ; 

(b) w(0) = w(l) = and w(u) < for all u G (0, 1); 

(c) w(u) > for allueX \ [0,1]. 

Remarks, (i) The proof of Lemma [7] in [KRllj is based on the analysis of the qualitative properties 
of the generalized Riccati equations in §4§ . 

(ii) The function u i— > w(u) from Lemma [7J can be extended naturally to a lower semicontinuous 
function w : K — )■ (— oo, oo] by w (M \ I) = oo. Then the extension, again denoted by w(u), has the 
following properties: 

• w is convex with effective domain T> w = X and X = {u £ R : w(u) < oo}; 

• the maximality of I implies that for u £ M\X there exists now* £R such that R(u, w*) = 0. 

The next theorem, proved in [KR11, Theorem 3.4], describes further properties of the function 
u i — y w(u) and specifies its relationship to the limiting cumulant generating function h defined 
in ([23]). 

Theorem 8. Let (X, V) be a non- degenerate affine stochastic volatility process that satisfies as- 
sumption A5 and let w(u) be given by Lemma^ Then the function h{u) defined in (|23p satisfies 

(24) h(u) = F(u,w(u)) for any u £ J := {s G I : F(s, w(s)) < oo}, 

where F is defined in ([5]) (see also ©J. Furthermore the inclusions hold: [0,1] C J C X. The 
functions w(u) and h(u) can be extended uniquely to cumulant generating functions of infinitely 
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divisible random variables and 



(25a) 



lim ip(t, u, 0) = w(u) for all it G X; 



(25b) 



lim —(p(t,u,0) = h(u) for all u G J . 



Remark. Since w(u) and h{u) can be extended to cumulant generating functions of some (infinitely 
divisible) random variables it follows that: 

• w (resp. h) is continuously differentiate in the interior of I (resp. J); 

• either h"{u) > for all u G J° or h(u) = for all u G M (this follows from ([24j) , (b) in 
Lemma [7] and assumption A3) . 

We say that R explodes at the boundary if lim n _ > . 00 R(u n , w n ) = oo for any sequence ((u n , iu n )) ngN 
in the interior T)° R that tends to a point in the boundary of Dr (both the boundary and the interior 
of T>r are relative to the topology of M 2 ) or equivalently Vr is open. By Proposition [3] (A), the 
gradient VF = (d±F, &2F) is continuous on T>° F . Analogously to the one-dimensional case (see (c) in 
Definition [5]) , we say that F is steep if linin^oo ||VF(u n , u> n )|| = oo for any sequence ((u n ,w n )) neN 
in the interior T>° F that tends to a point in the boundary of T>p. It is clear that if F explodes at 
the boundary, it is also steep but the converse may not be true. 

Before we state and prove the main results of this section (Theorem [TO] and Corollary [TT|) , we 
establish Lemma [9] which states that in an affine stochastic volatility model, the limiting cumulant 
generating function h cannot be identically equal to zero. This property will play an important 
role in understanding the limiting behaviour of the implied volatility smile (see e.g. Theorem 1 13|) . 

Lemma 9. Let (X, V) be a non- degenerate affine stochastic volatility process that satisfies assump- 
tion A5 and let h be given by (|23p . Assume further that the interior T>° F of the effective domain of 
F contains the set {(0,0), (1,0)} and that F(u,w) ^ for some (u,v) G T>p- Then h{u) > for 
all u G J ' \ [0, 1] and h{u) < for all u G (0, 1). Furthermore we have h"{u) > for all u G J° . 

Proof. Note that since h can be extended to a cumulant generating function of a random variable 
by Theorem [HI it is smooth in J° . Since h is either identically equal to zero or strictly convex on 
J by the remark following Theorem [8j the statement h"(u) > for all u G J° follows if we prove 
that h{u) < for some u G (0, 1). 

The function u i— > F(u, 0) is convex by (B) of Proposition [3j Furthermore it is either (I) strictly 
convex or (II) identically equal to zero (by A3). We analyse both cases. 

(I) Strict convexity and A3 imply that for u G (0, 1) we have F(u, 0) < 0. The same argument 
implies that for u G M \ [0,1], such that (u, 0) G V F , the inequality F(u, 0) > holds. The 
Levy-Khintchine representation of F in ([6]) implies that 



(26) 
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for any point in the interior of the effective domain T>p. It is clear from (I26p that &2F > on T>° F . 
Lemma [7] implies that for u G (0, 1) we have w(u) < 0. Identity (|24p in Theorem [8] yields 



The last inequality follows from the strict convexity of u \— > F(u, 0). If u G J° \ [0,1], then 
(u, 0) G T>° F and an analogous argument implies that h(u) > 0. The inequality at the boundary 
points of the interval J follows from the convexity of h. 

(II) Assume now that u i— > -F(-u, 0) is identically equal to zero. For any (u, 0) in the interior of 
the effective domain of F, the Levy-Khintchine representation of F in ([6]) yields 



This implies an = and m(d£) = (5q <X> z^)(d£), where ^(d^) is a Levy measure on (0, oo) with 
integrable small jumps and <5o is the Dirac delta. The condition F(1,0) = in A3 and the repre- 
sentation of F in ([6]) yield b\ = 0. Hence we have 



-%o\{°} v 7 

Since by assumption there exists (n, v) G Dp such that -F(ii, if) 7^ 0, either 62 > or v 7^ holds. 
Therefore identity ()24|) in Theorem [8l Lemma [7] and this representation of F conclude the proof. □ 

Remark. The assumption {(0, 0), (1, 0)} C T>° F in Lemma [9] ensures that the interiors of the effective 
domains of F and h are non-empty. It may not be necessary for Lemma [9] to hold. However, the 
assumption is crucial in Theorem [TO] and hence does not restrict the applicability of Lemma [9] in 
our setting. 

Theorem 10. Let (X, V) be a non- degenerate affine stochastic volatility process that satisfies as- 
sumption A5 and suppose that the function w 1— > F(0, w), where F is defined in ([5]), is not identically 
equal to zero. If R explodes at the boundary (i.e. T>r is open), F is steep and {(0,0), (1,0)} C T>° F , 
then the function h{u) is well-defined by (|23p as an extended real number for every «eK and its 
effective domain is given by T>h = J ( see (IMD for the definition of interval J). Furthermore h is 
essentially smooth and the set {0, 1} is contained in the interior (relative toM) ofT>h- 

Corollary 11. Let (X, V) be a non- degenerate affine stochastic volatility process that satisfies 
assumption A5 and assume that w 1— > F(0,w) is not identically equal to zero. If either of the 
following conditions holds 

(i) fi has exponential moments of all orders, F is steep, andT)° F contains (0,0) and (1,0), 

(ii) (X,V) is a diffusion, 

then the function h is well-defined by (|23p for every u G M. with effective domain T>h = J ■ Moreover 
h is essentially smooth and {0, 1} C T>^. 
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Proof of Corollary \lll Note that either of the conditions ((!]) or (Jn]) implies that T>r = M 2 and 
hence T>r is open. Therefore (JTJ) and the assumptions of Corollary [11] imply the assumptions of 
Theorem 1101 If (Jn|) holds, then (X, V) is a diffusion and 

u 2 

F(u,w) = an— + b\u + 02W with 011,62 >0 and b\ G R. 

Clearly T>° F = M 2 contains the set {(0, 0), (1, 0)} and F is steep if 62 is non-zero. In the case 62 = 0, 
the map w 1— > F(0, w) is identically equal to zero, which contradicts the assumption in Corollary II li 
Thus Corollary [11] follows from Theorem [10] □ 

Proof of Theorem\l(A The proof of this theorem is in two steps. In step (I) we show that {0, 1} C J° 
and that, if we extend h\j by +00 to M. \ J , we obtain an essentially smooth convex function. In 
step (II) of the proof we show that the limit in definition (123j) exists for any u G M as an extended 
real number and that definition of h in (I23j) agrees for every m£| with the extension of h\j from 
the first part of the proof. 

Step (I). Throughout this step we abuse notation by using h to denote the extension of h\j to R 
described above. Theorem [8] and the remark following it imply that h is essentially smooth (see 
Definition [5]) if it is steep. We will prove the steepness of h at the right endpoint = sup{n : u G 
J7} of the interval J and show that 1 G J° . The left endpoint U- = infjii : u G S} and the fact 
G J° can be treated by a completely symmetrical argument. 

Let (ii n ) n( =N be a sequence in J° converging to u + . We use the shorthand notation w n = w(u n ) 
and = lim n _ >00 w n , where u 1— > w(u) is the function given in Lemma [7] (note that the limit w+ 
exists but may be infinite since u 1— > w(u) is a cumulant generating function of a random variable 
and J C X = V w ). Since u \-t w(u) is convex on V w , the value w + is independent of the choice of 
sequence (n n ) ngN . 

Claim 1. The inequalities u+ > 1 and > hold. 

Indeed, since i?(l,0) = by assumption A3, we get that (1,0) G T>r = V° R . Assume now that 
114. = 1. Then by Lemma [7] we have w+ = and (u+,w+) G T)° R . Since R is continuously 
differentiable in T>° R and c?2i?(l,0) = x(l) < by assumption A5, the implicit function theorem 
and Lemma [7] imply that u+ is in the set 2?^ = 1°. Since (1,0) G P^, there exists u G 1° such 
that u > u + and (n, 0) G T)° F . Identity (|24p in Theorem [8] therefore implies that h(u) < 00, which 
contradicts the definition of u + . Therefore it + > 1. Lemma [7J implies that the sequence (w n )neN is 
eventually (certainly when u n > 1) non-decreasing and strictly positive. This yields that w+ > 
and the claim follows. 

Discarding finitely many elements we may assume that u n > 1 and w n > for all n. If u + is 
infinite, it is not in the boundary of J and the steepness of h follows. If n+ is finite but is 
infinite, identity (f2T|) and the assumption that w \-t F(0,w) is non-zero imply linin^oo h(u n ) = 00. 
The steepness of h follows from the convexity of h. Therefore in the rest of the proof we can assume 

(27) u + G (l,oo) and w + G (0,oo) 
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without loss of generality. 

Claim 2. The following statements hold true: 

(a) if u G 1°, where I is defined in Lemma [3 then (u,w(u)) G V° R and 

(28) = dtR(u, w(u)) + d 2 R(u, w(u))w'(u); 

(b) if u G J° n (l,oo), where J is defined in TheoremjHJ then (u,w(u)) G V° F and 

(29) h'{u) = diF{u, w{u)) + d 2 F(u, w(u))w'(u). 

The statement in (jaj) follows from Lemma assumption Dr = T>° R and the chain rule. To prove 
the first statement in (jb|), note that u G J° D (l,oo) C X° and hence X° \ [0, 1] 7^ 0. Lemma 
therefore implies that the function w : J — > R is strictly convex with u;(0) = itf(l) = and 
therefore strictly increasing on J° D (l,oo). Pick u' G ^7° Pi (l,oo) such that u' > u and note that 
(u',w(u')) G 2?f (by the definition of J) and (it',0) G 2?f (by representation © and < w(u')). 
Assumption (1,0) G Vp in the theorem and the fact w(u) < w(v!) imply that the point (u, w(u)) lies 
in the interior of the triangle with vertices (u', w(u')), (1, 0), (u', 0) in the convex set Dp. Therefore 
(u,w(u)) G T> F . Equality ()29|) follows by the chain rule. This proves the claim. 
Claim 3. The following holds for any strictly increasing sequence (un^ntaN with limit u+'. 

(a) if u + = sup J = supX, then 

\w'(u n )\—>oo as n — > 00; 

(b) if u + = supj" < supX, then 

\\VF(u n , w n )\\ — > 00 as n — > 00. 

To prove the claim, assume that the conclusion of (jaj) does not hold. Since the sequence 
(w' (u n )) n <=n is non-decreasing by Lemma [TJ there exists a finite positive number, denoted by 
w/(ii-|_) > 0, such that lmiyj—^oo w'(u n ) = w'(u + ). Claim 2 (jaj), applied to u = u n , implies 
(u n ,w n ) G V° R for all n G N and hence by (j27j) (u+,u;+) is in the closure of Vr. However (u+,w+) 
cannot be in the boundary of T>r since R explodes at the boundary by assumption and it holds 
lim^oo R(u mWn) — (recall that R(u n , w n ) = for all n G N). Therefore (u + ,w + ) G V° R . The 
derivatives d±R,d2R are hence continuous at (u+,w+) and, in the limit as n — > 00, formula ([28]) 
and the fact w'(u+) > imply 

d 2 R{u + ,w + ) = — - — . 

Therefore either = d\R(u+, w+) = d 2 R(u + ,w + ) or both partial derivatives at (u+,w+) are non- 
zero. Suppose the former. For an arbitrary u G (0, 1), the convexity of R yields 

-R(u, 0) = R(u + , w+) - R(u, 0) < VR{u + ,w + ) ■ (u+ - u, w+)' = 0. 

Since R(u, 0) < (see assumptions A3 and A4), this leads to a contradiction. Hence d\R(u + ,w + ) 
and d 2 R(u + ,w+) are non-zero and related by the equality above. By the implicit function theorem 
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there exists an open interval iV containing u+ and a function w : N — > R, such that R(u, w(u)) = 
for all u G N. This contradicts the maximality of X and proves Claim 3(jaj). Note that under 
assumption of Claim 3(Jb]), (u+, u;+) must be a boundary point of T>f- Since F is steep this implies 
\\VF(u n , w n )\\ — > oo and the claim follows. 

Theorem [TU1 follows easily if \\VF(u n , w n )\\ — > oo as n — > oo. Indeed, assumption (f2T|) and 
Lemma [7] imply that the sequence w'(u n ) > is strictly increasing and positive for all large n. 
Since F(0,0) = F(1,0) = 0, Proposition (B) implies that <9iF(l,0) > 0. The Levy-Khintchine 
representation of F in ([6]) implies d2F(u,w) > for all (u,w) G T>p. Since the gradient of the 
convex function F is monotone on T>° F and (u n ,w n ), (1,0) G 2?^ for all n, we find 



If \d2F(u n , w n )\ — > oo as n — > oo, the steepness of h at ii+ follows from (|27p . (|30p and the fact that 
w'(u n ) is strictly positive and increasing. If |9iF(ti n , w n )\ — >• oo as n — > oo, then, since (?2-F > on 
P^, formula (I29p implies Theorem 1101 

If ||VF(u n , ty n )|| does not tend to infinity as n — >• oo, the following facts hold: (u + ,w + ) G 
(since (ti + , w+) is in the closure oiVp and F is steep), \w'(u n )\ — ^ oo as n — > oo (by Claim 3) and 
d2F(u,w) > for all (u,w) £ T>° F (by Levy-Khintchine representation ([6]) of F). The next claim 
plays a key role in the proof of steepness of h. 

Claim 4. If || X7F(u n , w n )\\ does not tend to infinity as n — > oo, then d2F(u + ,w + ) > 0. In 
particular there exists 5 > such that d2F(u n ,w n ) > 6 for all large n G N. 

Note first that d2-F(u+,u) + ) is well-defined since (u + ,w + ) £ V° F . If d2F(u + ,w + ) = 0, differen- 
tiation under the integral in ([6]) implies that 62 = and the support of m is contained in the set 
M x {0}. This would imply that F(0,w) = for all w G R, which contradicts the assumption in 
the theorem. Hence the claim follows. 

To conclude the proof of Step (I), it remains to note that equality (|29|) applied at u n together 
with Claim 4 yield the steepness of h in the case \\VF(u n , w n )\\ does not tend to infinity. 

Step (II). We now prove that for any u G the limit in (|23p is equal to +00. This will conclude 

the proof of Theorem [TUJ 

Let t n 4 and define h n {u) = j- log Ee uXtn for all u G R. We know that lim^oo h n (u) = h(u) 
for all u G J . Moreover by Step (I), h(u) is steep at the boundary of J and G J° . Since X t is 
infinitely divisible for all t > (see |DFS03t Theorem 2.15]), there exist random variables X n such 
that h n (u) = log Ee Uj ^" (i.e. /i n is the cumulant generating function of X n ). Therefore there exists 
a random variable X such that X n — > X in distribution and, if we define H{u) = log Ee"^ for all 
u G R, the equality H(u) = h(u) holds on J . Since H is a cumulant generating function, it is 



dtFiu^w^iun - l) + d 2 F{u n ,w n )w n > diF(l,0)K - 1) + 3 2 F(1, 0)w n > 0. 



Therefore by (|29p we obtain 



(30) 
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lower semicontinuous and convex, and in particular continuously differentiable in the interior of its 
effective domain T>n- But h is steep and hence non-differentiable at the boundary of J . Therefore 
it follows that T>jj = J and H (u) = oo for all u G R \ J . However for all u £ R\ J , the Skorokhod 
representation theorem and Fatou's lemma imply 

liminf Ee uXn > Ee uX = e H(u) = oo. 



n— >oo 



Hence the equality lim n ^oo j- log Ee uXtn = oo holds for u G R \ J ■ This concludes the proof of 
Theorem [101 □ 

4.2. Degenerate affine stochastic volatility models. The remark following Definition [2] implies 
that in the case of a degenerate affine stochastic volatility process (X, V), the model S = e x is an 
exponential Levy model (note also that A5 in this setting fails). Therefore Definition[2]and imply 
that the characteristic exponent h(u) := F(u, 0) of X possesses a Levy-Khintchine characteristic 
triplet (5, a 2 , u), where 5, a 6 R and v a Levy measure on R \ {0}, and satisfies 

(31) h(u) = log E [exp {uX{) \X = 0] 

= u5+ \a 2 u 2 + I (j* - 1 - u-^-] iz(dfc) 
2 Jm\{o} V ?f + 1/ 

for all u £ C where the expectation exists. The independence and stationarity of the increments 

of X imply that S is a martingale if and only if h(l) = 0, which is, in terms of the characteristic 

triplet (5,a 2 ,u), equivalent to Jm j*,) e Kd£i) < 00 an d 

(32) 5 = -\a 2 - f (V - 1 - J±-) Kd6). 

2 Jm\{o} V ?i + 1/ 



The limiting cumulant generating function for the family of random variables (Xt/t)t>i, defined 
by the limit in (I18p . is in the case when X is a Levy process given trivially by h in (|3ip . which 
therefore also coincides with definition (|23|) . The martingale condition for S = e x and the convexity 
of h imply that [0, 1] is contained in the effective domain T>^- In the case of affine stochastic 
volatility models we had to establish Theorem 1101 to obtain sufficient condition for the set {0, 1} to 
be contained in the interior of the effective domain of h. In the setting of Levy processes it is 
well known (see e.g. [Sat99L Theorem 25.17]) that {0, 1} C T>° h if and only if 

/ e U0? V(d£i)+ / e ui ^u (d£i) < oo for some u < and u x >l. 

J (-oo,-l) J{l,oo) 

Condition (|33|) implies that the interior of the effective domain of h is of the form T>° h = (u_,ii + ) 
for some U- £ [— oo,0) and u + G (l,oo]. It is therefore clear that h is steep if and only if 

(34) f |£i|e* lU -i/(dfi) = oo and f ^ lU +u (dft) = oo, 

J (-oo,-l) J(l,oo) 

where the integrals are taken to be infinite if the integrands take infinite value for some finite £1 
(e.g. if ii_ = —00 or u + = 00). Note also that under assumption ([32]) . the Levy process X is non- 
constant if and only if there is a Brownian component (i.e. a 2 > 0) or its paths are discontinuous 



(33) 
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(i.e. v ^ 0). Hence the equality 



h"(u)=a 2 +[ £?e< M d6), ueV o^ 
Jr\{o\ 



implies h"(u) > for all u G T>° h . These arguments therefore imply Proposition I12( which is the 
analogue of Theorem [TOl for Levy processes. 

Proposition 12. Let X be a non-constant Levy process (i.e. the first component of a degener- 
ate affine stochastic volatility process) with state-space M, characteristic triplet (5, o~ 2 ,v) and the 
characteristic exponent h given by (|31|) , Assume further that conditions (|32|) . (|33|) and (|34|) are 
satisfied. Then the interior of the effective domain of h is an interval (u^,u + ), where u_ < u + 
are extended real numbers, h is a convex essentially smooth limiting cumulant generating function 
for the family {Xt/t)t>\ and the set {0,1} is contained in the interior ofT>h- Furthermore, h is 
smooth on T>° h and h"(u) > for all u G T>^. 

5. Rate functions and the option prices far from maturity 

In this section we describe the limiting behaviour of a family of European options under an affine 
stochastic volatility model S = e x . These results will be used in Section [6] to prove the formulae 
for the limiting implied volatility smile. 

In order to understand the limits of the vanilla option prices far from maturity in an affine 
stochastic volatility model S, we will need to apply the large deviation principle for the family 
(Xt/t)t>i under a risk-neutral measure P and under the measure P, known as the share measure]^ 
Recall that for every t > the measure P is equivalent to P on the u-field T% and the Radon- 
Nikodym derivative is given by 

dP x 
— = e * 
dP T t ■ 

The limiting cumulant generating function for (Xt/t)t>o under P is defined by 
(35) h(u) := lim — log E \e 

t—toc t 



uXt 



The function h and its effective domain D-^ satisfy 

(36) h(u) = h(u + 1) and V~ h = {u G R : 1 + u G V h }, 

where h is the limiting cumulant generating function defined in (|23|) and T>h is its effective domain. 
Note that G T>~ if and only if 1 G D° h . The identity in (|36p implies the following relationship 
between the Fenchel-Legendre transforms (see (|19j) for the definition) of h and h: 

(37) h*(x) = h*{x) - x for all x G R. 

Theorem [13] below describes the limiting behaviour of certain European derivatives under an 
affine stochastic volatility process (X, V). Before we state it, we collect the following facts. 



^The name stems from the fact that under P the numeraire asset is the risky security S 
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Remarks, (i) If (X, V) is a non-degenerate affine stochastic volatility process that satisfies assump- 
tions of Lemma [9] then the limiting cumulant generating functions h and h (defined in fj23[) and (|35p 
respectively) are strictly convex with strictly positive second derivatives in the interior of their re- 
spective effective domains. Remark (B) after Definition [5] implies that their convex duals h* and 
h* are strictly convex and differentiable with respective unique global minima attained at 

(38) x* = ti(0) and x* = h'(0) = h'(l). 
Lemma also implies the following inequalities: 

(39) x* < < x*. 

(ii) If (X, V) is a degenerate affine stochastic volatility process that satisfies assumptions of Proposi- 
tion [T2l then h! is strictly increasing on T>° h and its image is equal to R. The unique global minima 
of the Fenchel-Legendre transforms h* and h* are (by Remark (B) after Definition [5]) explicitly 
given by 

(40a) x* = h'(0) = -\a 2 - [ (e* - 1 - £i) K^i), 

(40b) x* = h'(l) = \o 2 + / (e^ (6 - 1) + l) Kd£i). 

1 iR\{0i v 7 



Formulae (|40p show that the inequalities in (|39p hold also in the degenerate case. 

(iii) In the case of the Black-Scholes model (i.e. v = 0), the assumptions of Proposition 1121 are 

satisfied. The effective domains of /ibs an d /^bs are equal to R and the following formulae hold 

(41) h B s(u) = ^a 2 (u 2 - u) and h B s{u) = ^cr 2 (u 2 + u) for u € R, 

(42) ^ s [x- a 2 ) = — (x+^j and ^ s (x; a 2 ) = — ( x - °— \ for x £ R. 
Therefore we have x* = —a 2 /2 and x* = a 2 /2. 

Theorem 13. Let (X,V) be a non- degenerate (resp. degenerate) affine stochastic volatility process 
that satisfies the assumptions of Theorem{Wi or Corollaries [77] (i), [77] (ii) (resp. Proposition 
Then the family of random variables {Xt/t)t>\ satisfies the LDPs under the measures P and P with 
the respective good rate functions h* and h* , where h is given in (|23|) (resp. (|3ip ) and h in (|35|) . 
Fix iEK, let x*,x* be as in (|38p (resp. (|40p ) and denote S = e x and y + := max{0, y} for y £ R. 
(%) T/ie asymptotic behaviour of a put option with strike Soe xt is given by the following formula 

v +-ii cTcc xt q\+] j x-h*(x) ifx<x*, 

hm t logE (,Soe - 5 t ) = < 

L -I [ x if x > x*. 

(ii) The asymptotic behaviour of a call option, struck at Soe tx , is given by the formula 

J — h* (x) if x >x* , 



lim t 1 log E 

t— ¥oo 



(St - S e= 



ifx<x*. 



limt- 1 logf5o-E[( ( St- ( S e a:t ) + l) = < 
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(Hi) The asymptotic behaviour of a covered call option with payoff St — (St — Soe tx ) + is given by 

if x >x*, 

x — h* (x) if x £ [x* , x*] , 
x if x < x* . 

Furthermore the convergence in (i)-(iii) is uniform in x on compact subsets o/R. 

Remarks. (I) The formulae in (i), (ii) and (iii) of Theorem 1131 are continuous in x since the value of 
the Fenchel-Legendre transforms h* (resp. h*) at x* (resp. x*) is equal to zero. Note further that 
the formulae in Theorem 1131 are independent of the starting value (Xq, Vq) of the model. 
(II) The reason for studying the limiting behaviour of the put, call and covered call in Theorem [T3l 
lies in the fact that these payoffs yield non-trivial limits on complementary subintervals of R, thus 
obtaining a non-trivial limit for every x £ R. This limit will be compared in Section [6] with the 
corresponding limit in the Black-Scholes model, which will yield the formula for the limiting implied 
volatility smile under afhne stochastic volatility models. 

Proof. Assume first that (X, V) is a non-degenerate affine stochastic volatility process. The limiting 
cumulant generating function h satisfies ()18j) and is essentially smooth either by Theorem 1101 or by 
Corollaries 1 1 1 1 (i) . 1 1 1 1 (ii) . Hence its convex dual h* is non-negative (by (|19j) and the fact /i(0) = 0), 
has compact level sets (by (|20p and E X>^) and is differentiable on D^* = R with strictly increasing 
first derivative (by Lemma [9] and Remark (B) following Definition [5]) . Therefore by Theorem [6] the 
family (X t /t)t>\ satisfies the LDP under P with the good rate function h* . Since 1 G V^, by ([36]) 
the function h satisfies the condition in (|18p . Therefore all of the assumptions of Theorem [6] hold 
under P and hence (Xt/t)t>\ satisfies the LDP with the good rate function h*. Furthermore h* 
enjoys the same regularity on 2?^„ = R as the rate function h*. The LDPs in the degenerate case 
follow from the same argument with Theorem [TUJ Corollaries [TT] (i) , [TTJ (ii) and Lemma replaced 
by Proposition 1121 

We now prove the formulae in Theorem [131 Without loss of generality we may assume that 
So = 1, i.e. Xq = 0. The following inequality holds for all t > 1 and e > 0: 



c 



tx (1 - e-) I {Xt/t<x - £} < (e xt - e^) + < e tx I {Xt/t<x} . 



Hence by taking expectations, logarithms, multiplying by 1/t and applying the LDP for (Xt/t)t>\ 
under P we obtain the following inequalities 



1 



y<x—e t—>oo t 



x- inf h*(y) < liminf-logE (e xt - e Xt Y <limsup-logE (e xt - e Xt V < x - inf h*(y). 



1 



y<x 



t— >oo t 

Since h* is continuous on R, strictly decreasing for x < x* and takes value at x*, the formula in 
Theorem [13] (i) holds. 

We now consider the call option case. The following inequality holds for all £ > 1 and e > 0: 



c 



Xt '1 ~ e" £ ) I { x t/ t> x+ e} < (* Xt ~ e*T < 



{X t /t>x}- 



LARGE DEVIATIONS IN AFFINE STOCHASTIC VOLATILITY MODELS WITH JUMPS 



23 



Again by taking expectations, changing measure to P, applying logarithms, multiplying by 1/t and 
applying the LDP for (Xt/t)t>i under P we obtain the following inequalities 



inf h* (y) < lim inf — log E 

y>x+e t— >oo t 



(e x t-e xt )~ 



< lim sup — log E 

t— >oo t 



,-X* xt 



e xt y <- inf h*{y). 

y>x 



Note that x* is a global minimum for h* at which value is attained. The continuity of h* implies 
the formula in Theorem [T3l (ii). 

In the case of the covered call, the following simple inequalities hold for all x G R: 



(43) 


e 1 {x t /t>x) 


< 


e x t _ 


(e x * 


-e tx ) + 




e Xt 


(44) 


eXtI {X t /t<x} 


< 


e x t _ 


(e x < 


-e tx ) + 


< 


e x < 


(45) 


e 1 {X t /t>x} 


< 


e x t _ 


(e x * 


-e tx ) + 


< 


e xt . 



T -{x t /t<x} +e xt I{x t /t>x}, 



Inequality (|43p and the LDP under measures P and P imply the inequalities 



e^P [X t /t >x\ < 1 - E 



(e Xt -e xt ) 



P[X t /t < x]+e xt P[X t /t > x 



< exp ( -t inf h*(y) + et + e xt exp ( -t inf h*(y) + et 
V y< x J V y> x 

for any x G K, e > and t large enough. Assume now x G [x*,x*] and note that in this case we 
have mi y > x h*(y) = h*(x) and inf^^^ h*(y) = h*(x). By ([37]) we obtain 



x + - log P [X t /t > x) < — log ( I - E 



( e X t _ e x t] 



< x - h*(x) + e + j log 2 



for any e and all large t. Therefore we find the inequalities 



x — h*(x) < liminf — log ( 1 — E 



(l_ E [(e x '-e^) + ]) 



< lim sup - log ( 1 — E 

t— >oo t 



( e X t _ e x t) + 



< x — h*(x) + e 



for all e > 0. This proves the formula in Theorem 1131 (iii) for x G 

Assume that x > x* and take expectations, change measure to P, apply logarithms and multiply 
by 1/t the inequalities in (I44p to obtain the following: 



i log P [Xt/t <x] < ilog( i -E 



( e x t _ e x t] 



< 0. 



Since mi y<x h*(y) = the LDP for (Xt/t\>\ under P implies the formula in Theorem 1131 (iii) that 
corresponds to x > x* . 

Finally let x < x* . Inequalities in (j45|) imply the following 



x + - log P [X t /t > x] < - log ( 1 - E 



( e X t _ e x t] 



< X. 



An application of the LDP for (Xt/t)t>\ under P completes the proof of part (iii). 
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We now show that the limits in the theorem are uniform in x on compact sets in R. Since the 
argument is similar in all the cases, we concentrate on Theorem [T3l (i). Let (xo,yo) be a finite 
interval in R and define for any 16R and t > 1 

V{t,x) = t _1 logE 

where v(x) denotes the continuous limit that appears in Theorem 1131 (i). It follows that 

V(t,x ) + v(x ) - v(x) <V(t,x) <V(t,y ) + v(y ) - v(x) for any x £ (x ,y ). 
We therefore find 

(46) \V(t,x)\ < max{\V(t,y )\,\V(t,x )\} + max{\v(x)-v(x )\,\v(x)-v{y )\}. 

Since we have already proved that limt_ ) . 0O \ V(t, x)\ = for every x and the limiting function v(x) 
is continuous, and hence uniformly continuous on every compact set, the inequality in (|46p implies 
that the logarithms of the put option prices converge to v(x) uniformly in x on compact sets in K. 

□ 



'Xt 



v(x), 



6. Asymptotic behaviour of the implied volatility 

The value C(Sq, K,t,o~ 2 ) of the European call option with strike K and expiry t in a Black- 
Scholes model (i.e. degenerate affine stochastic volatility model without jumps, see Section H~2j) is 
given by the Black-Scholes formula 

(47) C(S ,K,t,a 2 ) = S N{d+)-KN(d^), where d± = lo ^ S °/ K ^ aH / 2 

and iV(-) is the standard normal cumulative distribution function. Let S = e x be an affine sto- 
chastic volatility model from Definition [2] with the starting point Sq = e x ° . The implied volatility 
in the model S = e x for the strike K > and maturity t > is the unique positive number a(K, t) 
that satisfies the following equation in the variable a: 

(48) C{S ,K,t,a 2 ) =E[(S t -K) + }. 

Implied volatility is well-defined since the function a \— > C (So, K,t,o~ 2 ) is strictly increasing for 
positive a (i.e. vega of a call option ^(Sq, K, t, a 2 ) = SoN'(d + )\/t is strictly positive) and the right- 
hand side of (|48p lies in the image of the Black-Scholes formula by a no-arbitrage argument. Put-call 
parity, which holds since S = e x is a true martingale, implies the identity P (So, K,t,a(K,t) 2 ) = 
E [(K — St) + ] , where P (So, K, t, a 2 ) denotes the price of the put option in the Black-Scholes model 
with volatility a. 

If the affine stochastic volatility process (X, V) satisfies the assumptions of Theorem [T3"| then 
the implied volatility has the following limit 



(49) lim d(K,t) = 2^2/i*(0) = 2y/-2h ((/i')" 1 ^)) 

t— >oo 
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for any fixed strike K > 0, where h* is the rate function of the model (the second equality in (149 p 
follows from (|2ip and (|22|) ). Tehranchi [Teh09j proved that the first equality in (|49p holds uniformly 
in K on compact sets in R^o for non-negative local martingales with cumulant generating functions 
that satisfy certain additional conditions. Note that the limit in (|49p is independent of K, which 
corresponds to the well-known flattening of the implied volatility smile at large maturities. The 
uniform limit (in K) on compact subsets of R^Oj given in (|49p . also follows from Theorem 1141 for 
affine stochastic volatility processes (both in non-degenerate and degenerate, i.e. Levy, cases). 

In order to obtain a non-trivial limit at infinity we define the implied volatility at(x) for the 
strike K = Sq exp(tx), where x G R, by 

(50) (Jt{x) = o (So exp(tx), t) . 

We will show that if (X, V) satisfies the assumptions of Theorem [131 then the limiting implied 
volatility takes the form 



(51) CToo (x) = V2 



sgn(x* — x)J h*(x) + sgn(x — x*)y/h* (x) 



for x S 



where h* and h* are the Fenchel-Legendre transforms (see (|19p for definition) of the limit- 
ing cumulant generating functions h and h of X under P and P respectively and x* = h'(0), 
x* = h'(0) = h'(l). The function sgn(x) is by definition equal to 1 if x > and —1 otherwise. 

Remarks, (i) Under the assumptions of Theoremll3( the points x* , x* are the locations of the unique 
global minima of the good rate functions h* and h* respectively and by (|39p satisfy x* < < x* . 
Note that h*(x) < h*(x) for x > and h*(x) > h*(x) for x < and hence the following strict 
inequality £Too(x) > holds for all x E R. 

(ii) The function <7oo : M. — > (0, oo) in (|5ip is chosen so that it satisfies the following identities: 

(52) h% s (x;aoo(x) 2 ) = h*(x) and h^ s (x;a OQ (x) 2 ) = h*(x), 

where the polynomials /ig S (x;cr 2 ) and /ig S (x;<7 2 ) are given in (|42p . Quantities of importance in 
the proof of Theorem [J3] will be the following partial derivatives 

(53) ^ (x; a 2 ) = ^ (x; a 2 ) = ±, (a* + 2x) [o> - 2x) . 



(iii) In formula (f5T|) . the function o"oo(x) is given as a linear combination of \J h*(x) and y h*(x). 
The coefficients in this linear combination are not uniquely determined by identities (|52p (there are 
four possibilities). However definition (|5ip is the only choice that implies the following important 
properties 

(54) <7oo(x) 2 < 2\x\, for x G R \ [x* , x*] and <Joc{x*) 2 = 2x* , (Too(x*) 2 = 25?*, 

(55) (Too(x) 2 > 2\x\, for xG(x*,x*), 
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which will be crucial in the proof of Theorem [TH Note that (154[) and (I55p trivially hold in the 
Black-Scholes model: x* = —a 2 /2, x* = o~ 2 /2 and the limiting smile is constant and equal to 
a. The inequality in (|54p on the interval (— oo,x*) follows from the identity (Joo(x) 2 + 2x = 
4 h*(x) — ^Jh*(x) 2 — xh*(x) and the fact that for x < x* we have h*(x) > 0. Likewise the 

identity a^x) 2 + 2x = 4 h*(x) + ^Jh*{x) 2 - xh*(x) for x G (x* , 0) yields half of ([53]). The other 
half of (|54p and (|55|) follow from analogous identities involving /i* . 

(iv) The arbitrary choice for sgn(0) = 1 in (|5ip is of no consequence since h*(x*) = h*(x*) = 0. 

Theorem 14. Let (X, V) be an affine stochastic volatility process that satisfies the assumptions of 
Theorem{W[ Then we have 

lim at(x) = <7oo(x) for any 

t /*oo 

where o~t(x), defined in (|50p . is the implied volatility in the model S = e x for the strike K = Soe xt , 
and o~oo(x) is given in (|5ip . Furthermore for any compact subset C inM \ {x*,x*}, where x*,x* 
are defined in (|38p . we have 



sup 



o-t(x) 



as t /* oo. 



Remark. Theorem 1141 implies formula ()49j) obtained in [Teh09] : define x = log(K/So)/t and apply 
the uniform convergence of o~t(x) on a compact neighbourhood of zero. 

Proof. We can assume without loss of generality that Sq = 1. Assume further that xq > x* and 
pick e > 0. The goal is to find 5 > such that the following inequality holds for all large t: 

(56) \o~t(x) 2 — a oa {x) 2 \ < e, where x G (xq — 5,xq + 5). 

Inequality (}54"P implies that there exists e' G (0, e) such that <T 00 (xo) 2 +e / < 2xo and < o- 00 {xQ) 2 —e'. 
By ([53]) we conclude that a 2 h-> /ig S (x; a 2 ) is strictly decreasing on the interval (0, 2x) and hence 
obtain the following inequalities: 

{x^Oooixo) 2 - e') > h% s (xo^ooixo) 2 ) 



h 



BS 



> 



/ig S (x ; aoo(x ) 2 + e') . 



Since all the functions are continuous and identitiy (|52p holds, there exists a 5 > such that 
xq — 5 > x* and the strict inequalities hold 



(57) 



- /ibs (x;ctoo(x) 2 - e') < -h*(x) < -h% s (x; aoo{x) 2 + e') 



for all x G (xo— 5, xq+5). Theorem[T3]implies that the call option converges uniformly on the interval 
(xo — 5,xq + 5) to —h*(x) = limj^oo t^ 1 log E (e Xt — e xt ) + . In particular in the Black-Scholes 
model we get — /ig S (x;cr 00 (x) 2 ± e') = lim^oo t~ 1 log C(l, e x *, t, o"oo(x) 2 ± e') and the convergence 
is uniform in x on (xo — 5, xo + 5). Since o~t(x) satisfies E (e Xt — e xt ^ + = C(l,e xt ,t,a t (x) 2 ) by 
definition, the inequalities in ([57|) imply that 

C(l,e et ,t,v 0O (x)*-e!) < C(l,e xt ,t,a t (x) 2 ) < C(l, e xt , t, a^x) 2 + e') 
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for all x £ (xq — 5, xq + 5) and all large t. Since the Black-Scholes formula is strictly increasing in 
a 2 (i.e. vega is strictly positive), these inequalities imply (|56p . This proves uniform convergence on 
any compact subset C of (x*,oo). The proof for a compact set C C (— oo,x*) \ {x*} is analogous. 

We now consider convergence at the point x*. Pick any e > such that 0"oo(x*) 2 = 2x* > e. 
Identity (jUJ) implies that 

^BS {x*;voo(x*) 2 -e) > ^bs (x*;o-oo(x*) 2 ) = h*(x*) = < /i BS (x*; <7oo(x*) 2 + e) . 

The first inequality and the argument above imply that <7oo(5?*) 2 — e < <7t(x*) 2 for all large i. Since 
o~oc(%*) 2 + £ > 2x*, the second inequality, Theorem [131 yields 

1-C(l,e srt ,t,a t (x*f) > l-C{l,e**\t,a^{x*) 2 + s) 

for all large t. This implies <7j(x*) 2 < a^x*) 2 + e and hence proves the theorem for x*. The point 
x* can be dealt with analogously. □ 

The following corollary is a simple consequence of our results. 

Corollary 15. Let (X, V) be a non- degenerate affine stochastic volatility process that satisfies the 
assumptions of Theorem \13l Then there exists a Levy process Y , which satisfies assumptions of 
Theorem as a degenerate affine stochastic volatility process, such that the limiting smiles of the 
models e x and e Y are identical. 

Proof Let h be the limiting cumulant generating function for (X,V). Theorem [8] implies that h 
is a cumulant generating function of an infinitely divisible random variable. By Theorem [TP] the 
characteristic triplet of h satisfies conditions (|33|) and (fM|) . Therefore, if we define a Levy process 
Y with this characteristic triplet, Theorem 1141 and formula (|5ip imply that models X and Y have 
identical limiting volatility smiles. □ 

Remarks, (i) In other words Corollary [15] states that in the limit, non-degenerate affine stochastic 
volatility models cannot generate the behaviour of implied volatility, which is different from that 
generated by the processes with constant volatility and stationary, independent increments, 
(ii) Corollary [15] suggests the following natural open question: can any limiting smile of an expo- 
nential Levy model be obtained as a limit of implied volatility smiles of a non-degenerate affine 
stochastic volatility process? It is not immediately clear how to approach this problem because 
the characterisation of the limiting cumulant generating function h of a model (X,V), given in 
Theorems [8] and [10] does not give an explicit form of Levy-Khintchine triplet of h. 

6.1. Examples of limiting smiles. We now apply the analysis to the examples of affine stochastic 
volatility models described in Section 12.21 In each of the cases the limiting cumulant generating 
function h is available in closed form. If the assumptions of Theorem QJJ] or Corollaries llll(i).llll(ii) 
are satisfied, then the convex dual h* is a good rate function and hence the formula in (|5ip defines 
the limiting smile as maturity tends to infinity. 
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6.1.1. Heston model. The characteristics F, R are given in (llip and x( u ) = u Cp ~ ^ ( see (H2p ). 
Assumption A5 is satisfied if and only if x(l) < 0) which is equivalent to A > Qp. Since \9 ^ 
it follows that w i— >■ .F(0, to) is not identically 0. Since the assumptions of Corollary [11] (h) are 
satisfied, Theorem 1141 implies that the limiting smile is given by the formula in (|5ip . where 

(58) h(u) = -—(x(u) + ^W)) and A(u) = X (u) 2 - CV - u). 

This implies the main result in [FJllj . [FJMllj . A first order asymptotic expansion for the large 
maturity smile in the Heston model was obtained in [FJM10] using saddle point methods. 

6.1.2. Heston model with state-independent jumps. The functions F, R are given in ()13p and x( u ) = 
u^p — A. As in Section [6.1.1[ A > (p implies that (A, V) defined Section [2.2.21 is a non-degenerate 
affine stochastic volatility model that satisfies A5. As before assumption X6 ^ implies that w h-> 
F(0,w) is non-zero. k(u), defined in (|14p . is a cumulant generating function of the compensated 
pure-jump Levy process J. Assume that there exists k_ < such that < oo for u > k_, 
|/5(u)| = oo for u < k_ and (|34|) holds for u_ = k_ and u + = oo (e.g. if the distribution of 
the absolute jump heights is exponential with parameter a > 0, then K- = —a). Under these 
assumptions on state-independent jumps, the function F in (|13ap is steep and {(0, 0), (1, 0)} C T>° F . 
Hence Theorem [10] implies that the limiting cumulant generating function is of the form 

X9 / \ 
Ku) = + V A (u)J +7t(u), 

where A(u) is as in (|58|) . and Theorem 1141 yields the limiting smile formula in (|51|) . Note also that 
condition (|34p on the jump measure is not necessary if A(n) < for some u > k_, since in this 
case F in (|13ap is automatically steep. 

6.1.3. A model of Bates with state- dependent jumps. The functions F, R are given in (|15p . Again we 
assume A > Qp and A# 7^ 0, which implies that (X, V) defined in Section [2.2.31 is a non-degenerate 
affine stochastic volatility model that satisfies A5. It is clear from (|15a|) that the assumptions of 
Theorem 1101 on F are satisfied. Let k(u) be as in (fl~5bj) and assume that either |k(u)| < 00 for 
all u £ R or there exists k_ G R such that < 00 for u > k_, = 00 for u < n- and 
lim^^oo = 00 for any sequence (u n )neN with « n | k_. Then Vr is open in R 2 (see (|15b|) ) 
and Theorem 1101 implies that the limiting cumulant generating function takes the form 

H u ) = -T2 (^( n ) + y /A ( u )) » where A ( n ) =x(u) 2 -C 2 (u 2 -u + 2k{u)). 

6.1.4. The Barndorff- Nielsen- Shephard model. The functions F, R are given in (117p , Since x( u ) = 
— A < and the jump measure is non-trivial (i.e. v 7^ 0), the process (X, V) defined in Section ^. 2, 41 
is a non-degenerate affine stochastic volatility process. Assume that k(u) defined in (|16p is either 
finite for all u £ R or there exists k+ > with < 00 for u < k+, |k(u)| = 00 for u> k + and 
and (|34p holds for n_ = k_ and u + = 00 (e.g. if the distribution of the absolute jump heights is 
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exponential with parameter a > 0, then k_ = —a). Then F (see (|17aj> ) satisfies the assumptions 
of Theorem [10] and the limiting cumulant generating function is of the form 

h(u) = Xk (^r+u(p - ^-Y) - u\k{ P ). 



\2X V 2A 

6.2. How close are the formula a oc (x) and the implied volatility at(x) for large maturity? 

In this section we plot the difference {cr^x) — cr t (x)\ for t S {10, 15} and x G [—0.1,0.1] for the 
models with jumps from Section 16.11 (see Figure [T]). In the case t equals 10 years the error is 
approximately 45 basis points (bp) with the strike K ranging from 30% to 200% of the spot. At 
the maturity of 15 years the error is approximately 20 bp and K ranges between 20% and 400% of 
the spot. 

In the cases of Heston with state-independent jumps and Bates with state-dependent jumps we 
took the following diffusion parameters 

A = 1.15, C = 0.2, = 0.04, p = -0A 

and the Levy measure v (d£i) = ae 01 ^ 1 I^ 1<0 jd^± with a = 0.6. The compensated cumulant gener- 
ating function k(u) (see ()14[) for the definition of k(u) in Section 12.2.21 and note that it takes the 
same form in Section I2.2.3[) is in this case given by 

~t \ u(u — i) . , , 

k(u) = - —. r for all u S (— a,oo . 

y ' (u + a)(a + l) v ' ' 

In the case of the BNS model we took a pure-jump subordinator J with Levy measure v (d^) = 
a6e _b ^ 2 /|^ 2> o}d^2- The cumulant generating function (|16|) is given by k(u) = au/(b — u) for u < b. 
We used the following values for the parameters 

a = 1.4338, b = 11.6641, A = 0.5783, p = -1.2606, 

which were taken from [Scli03t Section 7.3] where the model was calibrated to the options on the 
S&P 500. 
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Figure 1. This fi gure contains the plots of the function x i-> \o x (x) — <Jt(x)\ in the interval 
x G [—0.1, 0.1] for the models with jumps from Section [6T] and maturities t £ {10, 15}. The values 
of the model parameters used are given in Section [ 
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